Interaction between vortices in models with two order parameters 
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The interaction energy and force between widely separated strings is analyzed in a field theory 
having applications to superconducting cosmic strings, the SO (5) model of high-temperature super- 
conductivity, and solitons in nonlinear optics. The field theory has two order parameters, one of 
which is broken in the vacuum (giving rise to strings), the other of which is unbroken in the vacuum 
but which could nonetheless be broken in the core of the string. If this does occur, there is an effect 
on the energetics of widely separated strings. This effect is important if the length scale of this 
second order parameter is longer than that of the other fields in the problem. 
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In at least three very different contexts, nonlinear effects in field theories with spontaneous symmetry breaking give 
rise to topological solitons wherein a second order parameter (unbroken in the vacuum) attains a nonzero expectation 
value in the core of the soliton. 
| First, bosonic superconducting cosmic strings ,l| arise in a model with two U(l) symmetries. The first of these 

■ (which could be gauged or ungauged) is spontaneously broken by a complex scalar field tp; this gives rise to the 
possibility of string solutions, where the phase of ip changes by 2ir around a large loop in space. In the core of such a 
string, (p — > 0. The second symmetry is gauged and unbroken, and is identified with electromagnetism. It is supposed 
that an electromagnetically charged scalar field a exists; although its VEV is zero, the potential V(a, ip) is such that 
if one forces tp to zero, then V(a, 0) is minimized for a ^ 0. Thus, it is possible that a attains an expectation value 

• inside the core of a string, making the string a superconducting wire. Whether or not this actually occurs is a detailed 
dynamical question. A variety of generalizations of this idea have been discussed in Refs. @-G9- 

Second, in nonlinear optics, one can consider beams for which the amplitude of the electromagnetic field envelope 
function vanishes and its phase changes by an integer multiple of 2tt at a given point. In a nonlinear medium such 
a configuration can be a stable solution of the equations of motion, known as an optical vortex soliton (see [TH ] for 
references and review). By coupling a second propagating mode to the first, the vortex can act as a waveguide for the 
second mode, which is confined to the core of the vortex (for theoretical work, see 0, I f° r recent experimental 
results see |lj,[l5|,[l(j)- The first and second modes are analogous to the fields ip and a above, respectively, although 
.,— I | in the optical context both fields are ungauged. (It should also be mentioned that the fields are not really order 

■ parameters here.) 

Finally, the SO (5) model of high-temperature superconductivity (HTSC) 01 is a "unification" of the two phenomena 
that occur in these materials at low temperature: antiferromagnetism (AF) at low doping (including the undoped 
case) and superconductivity (SC) at higher doping. The Ginzburg-Landau (GL) model is written in terms of a 
five-component real order parameter composed of a three-component field rj describing antiferromagnetism and a 
two-component field <f) describing superconductivity. In the SC phase, the latter field attains an expectation value, 
and strings (vortices) exist, just as in conventional (SO(2)) superconductivity. It is possible that r\ attains a nonzero 
expectation value inside the core, making the core AF, which could provide an experimental test of the SO (5) model 
.18]- 

In this paper, we will study a field theory which applies (with minor modifications) to any of the above situations, 
and will study the energetics of two (widely-separated) vortices. In the cosmic string context this has an influence on 
the dynamics of a network of strings [T^. I20L uA\ ; in optics this affects the stability of solitonic waveguides [Tlj ; and 
in the superconducting case this is one way to study the type of superconductor described by a given model, which 



provides another test of the SO (5) model 
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In the next section, we will establish notation and review previous work (expressed in the language of superconduc- 
tivity, but easily translated into the other contexts), wherein the effect of the AF phase on the magnetic behaviour 
was studied, and the conditions under which SO(5) superconductors are type I or II were determined. This was done 
in two ways. First, the free energy F of a vortex as a function of its winding number m was calculated numerically 
|22|. If the energy per unit winding number T{m) = F(m)/m is a decreasing function of m, type I superconductivity 
results, since it is then energetically preferable for a given flux to penetrate a superconductor in one large region, inside 
which superconductivity is destroyed. On the other hand, if !F{m) is an increasing function of m, the superconductor 
is type II, since the energetics then prefers a given amount of flux to be divided into a network of vortices of winding 
number m = 1. Second, the type of superconductor was deduced from a determination |23| of the critical magnetic 
fields (see, for example, 24] ). While this approach is perhaps less intuitive than the previous one, it has the advantage 
that a simple analytical determination of the boundary between type I and type II superconductivity can be made. 

A third way of differentiating between type I and type II superconductors is through the force between two widely- 
separated vortices: they attract or repel for superconductors of type I or II, respectively. This approach is very 
intuitive, and is also applicable to both the cosmic strings and nonlinear optics contexts (where the notion of critical 
magnetic fields which restore the symmetry broken in bulk doesn't apply). In this article, we will examine this third 
approach, adopting for the most part the language of SO(5) superconductivity, though much of the discussion can 
easily be exported to other contexts such as those mentioned above. 

The method used is essentially that used by Speight |2^| . We write an expression for the energy of two vortices 
using a point vortex approximation, wherein the full nonlinear free energy is replaced by a linearization of it plus 
point sources for each vortex. The approximation is expected to be valid if the intervortex separation is much larger 
than the core size. The approximate free energy can be written as the sum of the individual vortex energies and an 
interaction energy, from which the intervortex force can easily be found. In conventional SC, vortices attract or repel 
one another in type I and II SC, respectively. As we will see, this is not always the case here: for certain values of the 
parameters of the SO(5) model, a superconductor which is type I (in the thermodynamic sense) will have repulsive, 
not attractive, vortices. This unusual behaviour is unlikely to be seen in HTSC, however, as it occurs only for small 
values of the GL parameter (k ~ o(l)), whereas in all known HTSCs, k > 50. 

II. REVIEW OF PREVIOUS WORK 

The model we wish to consider is described by the following two-dimensional free energy: 

+ ^(Vr,) 2 + n<^)}- (1) 

Here 4> and t) are the SC and AF order parameters, respectively. The former is a complex field associated with the 
U(l) gauge field A, while the latter is a real triplet whose SO(3) symmetry is ungauged. ^ 
The potential is taken to be an even, quartic function of \(f>\ and \f)\: 

v&t» = - a -^\4>\ 2 a \u 2 + \ (h 4 + m 2 H 2 + \v\ 4 ) ■ 

For simplicity, we have imposed SO(5) symmetry on the quartic couplings; the quadratic couplings are negative to 
give rise to symmetry breaking. The ground state depends on the value of the parameter (3 = a\/a\, and is SC if 
(3 < 1. We are primarily interested in this case; the ground state can be written (\<t>\, \r)\) = (v, 0), where v = ai/Vb. 

The model can be simplified somewhat by rescaling the fields and the position variable. Defining dimensionless 
(unhatted) quantities 

2 a\c\Jm* •> , / m*c 2 - a?c 2 m* 

we have the dimensionless free energy 

F = \ J d2s { (v x A)2 + h ( IHV - kA)0 ' 2 + (Vr?)2 ) 

_^_^ + i (02+r?2)2+ i| ; (2) 



F = t#A£±W + K 
J I 8tt 2rr 
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where 



r\ = |tj|, the derivatives are now with respect to s, and 



m c 




is the usual GL parameter. (A constant has been added to F so that the ground state has zero energy.) 

From (J2J, we see that the behaviour of the model is completely determined by two dimensionless parameters: k 
and the parameter /3, which is the ratio of the quadratic coefficients of 77 and <p. 

An observation which will be useful below is that if we set 77 = 0, then the free energy (a function of <j) and A only) 
reduces to that of the SO (2) model. Thus, for example, the energy of a static configuration for which 77 = in the 
SO (5) model will be exactly equal to the energy of the same field configuration in the SO (2) model with the same 
value of K. 

Since the complex field <p attains a nonzero VEV, vortex solutions exist, wherein the phase of cj> changes by 2-7T 
around a large circle in space. One can also consider configurations of higher winding number; a rotationally symmetric 
ansatz of winding number m ( "m- vortex" ) is essentially the conventional one for <f> and A together with a rotationally 
invariant ansatz for rj, whose orientation is taken to be a fixed, arbitrary unit vector e: 



0( S ) = <p{ a y 



Ms) 



'-A(s), t)(s) =&n(s). 



(3) 



As in the conventional case, these solutions carry a magnetic flux proportional to m (specifically, the dimensionless 
flux is $ = —2irm/K). 

As mentioned in the Introduction, a non-vanishing 77 field may appear inside the core of the vortex due to the 
competition between the potential and gradient energy terms. By examining the potential, one can see qualitatively 
that for fixed K, as (3 increases, the impetus for 77 to be nonzero in the core of a vortex grows, since the potential 
energy savings gained by having 7/ 7^ increases with (3, while the kinetic energy cost is independent of (3. Thus, we 
can define a critical value of f3, above which the core of a vortex is AF and below which it is normal. This critical 
value depends upon k, and also upon the winding number m of the vortex. (The dependence on m can also be argued 
qualitatively, by noting that the larger the winding number, the greater the region in which ip is nearly zero, and the 
greater the impetus for 77 to be nonzero.) (3af{k, m) was found numerically for m — 1 to 5 in Ref. |22l |. 

The fact that the model reduces to the SO (2) model if 77 = enables us to easily make contact with that model; we 
need only set (3 to zero, since then there is never any reason for 77 to be nonzero in all situations we will consider here. 

Substituting the ansatz © into (J2J, we find the following free energy: 



F(m) 




^ 2 +(^+^)V+7/ 2 



V-/V+o(¥> 2 +?7 2 ) 2 + 



(4) 



The free energy, written as a function of m, depends also, of course, on the two parameters of the model, k and (3. 
The equations of motion that follow from ||3} are 



1 



-ip' - [ — +kA) ip 



+ p(l - ^ 2 - T? 2 ) = 0, 



(5) 



(6) 



A" + -A'-^A- [j^ + A ^ = 0. (7) 
These equations cannot be solved analytically; however, asymptotic solutions for large s can be found. Defining 



(p( 8 ) = 1 - f( s ) and A(s) 



a(s), 



(8) 



the fields /, a and 77 approach zero exponentially as s — » 00. The linearized equations in these fields have solutions 
f(s)=C f K (V2Ks), a(s) = CaK^s), r)(s) = C v K (^/l - 0ks) (9) 
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FIG. 1: Stability of vortex lattice in the SO(5) model. Curves A and B separate the normal core region (left) from the AF 
core region (right) in the case of m = 1 and m — 2 vortices, respectively. Curve C depicts the border between type I and type 
II behaviour (see llll l. These curves divide the /3-k plane into sevelar regions, the most interesting of which are regions 5 and 
6, where the SC is type I yet vortices (of winding number 1) repel one another. 

where Kq and K\ are modified Bessel functions of the second kind and the C's are constants not determined by the 
linear equations. These results will be useful for deriving the potential energy between a pair of widely-separated 
vortices in the next section. 

By studying either the energetics of vortices as a function of their winding number |22j | or by studying the critical 
magnetic fields p3| . one can determine a curve in the /3-k plane indicating the boundary between type I and type II 
behaviour, which yields a surprising result. 

To see this, let us first recall the situation in the conventional SO (2) model, focusing our attention on vortex 
energetics. There, the only parameter is k, and k = l/v2 = n c is the critical value separating type I [k < k c ) and 
type II (n > k c ) behaviour. For example, one finds that the energy per unit magnetic flux of a vortex T(m) — F(m) / m 
decreases or increases with m according to whether k < k c or k > K c . In either case, T{m) is a monotonic function 
of m, tending towards a constant value for large m. 

In the SO (5) case, the monotonicity of J-(m) is no longer guaranteed. The reason is that for fixed k and /3, as m 
increases one can go from a normal vortex core to an AF core. As long as the core is normal, its energy is exactly 
as in the SO (2) model, as can be seen from @. However, as m increases, eventually the core becomes AF, at which 
point the energetic picture changes: in particular, the energy of all subsequent m-vortices will be reduced (relative 
to the equivalent SO(2) model), lessening the degree to which the superconductor is type II, and in some cases even 
changing the superconductor from type II to type I. For fixed k, this effect is stronger the larger the value of f3; we can 
therefore define a critical value (3 c (k) as that value for which T{m) goes from an increasing (type II) to a decreasing 
(type I) function of m, in the limit of large m. 

It is strange, however, that (in contrast to conventional SC) !F{m) is not necessarily a monotonic function: it could 
increase with m while the core is normal and decrease with m subsequently. The large-™ behaviour indicates a type I 
superconductor, yet the energy of a 2-vortex is more than double that of a single vortex, which is a feature of type II 
superconductivity. This does not cast doubt on the fact that such a SC is indeed type I, in a thermodynamic sense; 
however, it will certainly affect vortex dynamics. 

The behaviour at large m is confirmed by studying the critical ma gnet ic fields H c and H C 2- As in conventional SC, 
the relative magnitude of these is an indicator of the type of SC. In [23] the critical fields were found to be given by 
the following expressions: 

H c (/3) = ff c Vl-/3 2 , H c2 (P) = V2 K H° C (1 - 0), (10) 

where H® is the thermodynamic critical field in the conventional SC model. The border between type I and type II 
will then be given by equating these two fields, leading to the following expressions for /? c (k) and, by inversion, k c ((3) 
(that value of K above which a SC is type II, for fixed (3): 

&(*) = |£^I> «c(/?) = (1/^2)^/(1 + /?)/(! -/?)■ (11) 

This line, as well as the lines separating the normal/AF core boundary for m = 1 and m = 2 vortices, lines C, A, and 
B respectively, are shown in Figure ^ 
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III. INTERVORTEX FORCE 

As mentioned in the Introduction, a third way of determining the magnetic nature of a superconductor is through 
the interaction energy and force between vortices. We will study the interaction of widely-separated vortices in the 
SO (5) model. In conventional superconductivity, vortices of any separation are attractive or repulsive in the case of 
type I or type II superconductors, respectively, as was shown in great detail numerically by Jacobs and Rebbi |26| . 
We would thus expect the intervortex force in the SO (5) model to change from repulsive to attractive as we cross 
curve C in Figure ^ from left to right. This turns out to be not always the case. The reason is that vortex energetics 
does not always reflect the type of superconductor in the SO(5) model, as a comparison of the following two facts 
demonstrates. On the one hand, the type of an SO (5) superconductor depends critically on the AF sector (see |23jl. 
On the other hand, when vortices have normal cores, their behaviour (their interaction energy, in particular) is blind 
to the AF sector of the model, as was noted in the previous section. Thus, vortices with normal cores cannot be 
expected to necessarily behave in the fashion dictated by the type of superconductor involved. In particular, it is 
possible that vortices in a type I (in the sense of bulk thermodynamic properties) superconductor repel one another. 

One subtlety which arises in the interaction of vortices once their cores become AF is that while the orientation 
of rj can be taken fixed in the core of any given vortex, in the absence of anisotropics, there is no reason to expect 
the AF cores of different vortices to be oriented in the same direction. When it becomes energetically favourable to 
develop an AF core, each vortex randomly selects a drection in which rj will point. The interaction energy between 
vortices will depend on the angle between the two AF order parameters. As we will see, the contribution of r\ to the 
vortex interaction energy is most important if the two orders parameters are parallel. Our discussion will assume for 
the most part that this is the case. 

To clarify the situation, consider the interaction of ordinary vortices (that is, vortices of winding number 1). The 
core of these vortices is normal or AF to the left or right of curve A of Figure respectively. Their interaction energy 
will be shown to obey the following behaviour. If the core is normal (to the left of curve A), the energetics is exactly 
as dictated by the SO(2) model: attractive for n < 1/V2 and repulsive for n > l/\/2. If the core is AF (to the right of 
curve A), then (except for a very thin transition zone wherein the AF core turns on) the energetics correctly reflects 
the type of superconductor involved (repulsive or attractive to the left or right of curve C). 

If we imagine fixing k and increasing /3, then there are three cases. First, if n < 1/V2 nothing dramatic happens: 
for all (3 the model describes type I superconductivity and vortices attract. Eventually curve A is crossed, so the 
vortices develop AF cores; they remain attractive, however. 

Second, if k lies above its value at the point of intersection of curves A and C (k > 2.25), the situation is slightly 
more interesting: the core turns AF to the right of curve A, but the vortices remain repulsive until curve C is crossed, 
after which the superconductor is type I and the vortices attract. 

The third, intermediate case, where curve C is crossed before curve A as (3 increases (1/V2 < n < 2.25), is the most 
interesting. When curve C is crossed, the superconductor becomes type I, so the vortices ought to attract. However, 
their cores are normal, their energetics is as in the SO(2) model, and they repel, since k > 1/V2. This anomalous 
behaviour persists until curve A is crossed, when the core becomes AF and the vortices start to attract one another 
- the expected behaviour for a type I superconductor. 

In summary, curve C delineates the boundary between type I (to the right) and type II (to the left) superconduc- 
tivity; however, in regions 5 and 6 the superconductor has a sort of identity crisis: it is type I, but its vortices behave 
as type II vortices. |30j 

One can also study the interaction energy of vortices of higher winding number. The behaviour is qualitatively 
similar to that for 1-vortices, although the region of disagreement between vortex energetics and the type of super- 
conductor is reduced. For 2-vortices, for example, the core turns AF at a lower value of /3 for any given k (curve 
B rather than A), and consequently it is only in region 6 that one finds disagreement. This trend continues as the 
winding number is increased. 

With this summary of the main results concluded, let us describe the calculations involved in determining the 
interaction energy of a pair of vortices. For arbitrary separation, an analytic solution is not possible; even numerically, 
the problem is extremely difficult. However, the interaction between widely separated vortices is considerably more 
tractable, and is clearly of interest, in that it is an important ingredient in determining the dynamics of widely 
separated vortices, and is normally directly related to the type of superconductor involved. 

If the vortices are widely separated, one can argue that the details of the fields in the core ought to be unimportant, 
and we can model them by a simplifying approximation, following Speight 25]: we can consider a linearized theory 
with point sources added at the location of the vortices in such a way that the long-range fields produced by the point 
sources in the linear theory agree with those of the vortices in the original theory. 

To derive the linearized theory, we must expand the free energy around the asymptotic values of the fields. 
Before doing this, however, it is useful to eliminate the phase of <j), as this degree of freedom becomes the longitudinal 
gauge field. Thus, we take (f> = 1 — f real, and expand @ in powers of /, A and ry, up to quadratic terms. This 
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results in: 

Ffrcc = \ Jd 2 s |(V x A) 2 + A 2 + -1 (V/) 2 + 2f + 1 (Vry) 2 + (1 - /?)r, 2 | . (12) 

To this, we must add couplings to sources: 

Source - J d 2 s(pf + a ■ q + j ■ A). (13) 

The equations of motion which follow from -Ff rcc + -F S ource are 

V 2 / - 2k 2 f = K 2 p, (14) 

V 2 r/ - (1 - P)k 2 t] = K 2 cr, (15) 

V 2 A-A=j. (16) 

The sources are to be chosen to give rise to fields which coincide with the asymptotic vortex fields, after elimination 
of the phase of cf> via a gauge transformation. These asymptotic fields are (sec ©): 



f(s) = C f K {V2Ks), A i (s)=e ij ^C a Kx(s), rj(s) = eC v K (^l - 0ks). (17) 



We require sources such that the solutions of (|14lltj|) are l|17|) . For / and rj the answer is found directly by 
substitution once we recognize that 

(V 2 - fi 2 )K (ps) = -2irS 2 (s); (18) 

we find 

p(s) = -2^X<S 2 (s), o.( s ) = _ 2 ^e%J 2 (s). (19) 

K K 

For the gauge field, differentiation of (|18|l and substitution yields 

Ms) = 2iiC a e l3 d. J 5 2 (s). (20) 

The vortex is now described by point sources of magnitudes 2wCf/K 2 and 2-kC^/ 'k 2 for / and rj, respectively, and an 
infinitesimal current loop of magnetic dipole moment 2irC a for A. 

We are now in a position to derive the interaction energy of two widely separated vortices. Suppose the vortex 
positions are Si and S2. In the point vortex approximation, each vortex is described by sources of the form I|19l20fl . 
displaced to the position of the vortex. Linearity of the equations of motion then implies that the fields will just be 
the sum of the individual vortex fields, and the energy of the configuration will be given by -Ffree + -FLourcc, with fields 
(f,T), A) = (/i +f2,Vi + r l2i A i + A 2) and sources (p, tr,j) = (pi + p 2 ,cr 1 +er 2 ,ji +J2), where subscript 1,2 indicates 
s — » s — Si^. We can subtract off the energy of each vortex to obtain the following interaction energy: 

*i«t = / d 2 s ( Pl f 2 + a x ■ r, 2 + j! • A 2 ) = F f + F v + F A . (21) 

To evaluate Ff we simply substitute p\ = p(s — Si) and f% = f(s — S2); we find 

/( C \ G 2 

d 2 s f-27r-|5 2 (s-si)J C/Xo(V2k|s-s 2 |) = -2iT-^Ko(V2Kd), 

where d = \s% — s 2 \ is the separation of the vortices. The same argument applies to F v , yielding 



(22) 



C 



2 



F v = -27r-£ei ■ e 2 K (y/l - 0nd). (23) 
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FIG. 2: Behaviour of the potential in the case of the S0(2) model (/3 = 0) for various values of k: k = 0.6 (type I), n = n c 
(critical), and k — 1.0 (type II). Note that as d approaches zero, the exponential behaviour is lost and the curves can not be 
trusted anymore. 



For F a , a similar (but slightly more complicated) procedure leads to 

F a = f d 2 S (27rC a e lJ d S] 6 2 (s~s 1 ))eJ-^—^C a K 1 (\ S ~s 2 \) 

J I s ~ s 2| 

OSlj \ I Si - S2| / 

= 2nC 2 a K (d). (24) 

Thus, the interaction energy of two widely separated vortices takes the form 

/ Q2 qi \ 

U{d) = 2tt C 2 a K {d) - -tK {V2nd) - -|ei • e 2 K (y/l - p K d) . (25) 

The first two terms give the interaction energy in the SO(2) case, and coincide with the results of |27|. The 
third term is the effect of the AF cores. As can be seen, it depends on the angle between the two order parameters. If 
they are parallel the contribution to the interaction energy is maximally negative, resulting in the greatest attraction, 
while antiparallel order parameters yield a positive (repulsive) force. If the order parameters are orthogonal the AF 
contribution to the interaction energy is zero and the results coincide with the SO (2) case. Since the parallel case is 
the most interesting and is also energetically preferred, we will assume in what follows that the order parameters are 
parallel, that is, that &i ■ e 2 = 1. 

In order to obtain useful information from 125|) . we must determine the constants C a , Cf and C„ by comparing @ 
with numerical solutions of the nonlinear equations 1)5171) , which have been obtained previously [221 123| . 

We will first study the conventional (SO(2)) case to make contact with previous work; This can be achieved within 
the SO (5) model by taking f3 = 0. Subsequently, we will examine the general SO (5) case given by 1)25(1 . The coefficients 
C a , Cf are displayed for a variety of n (and for j3 = 0) in Table 1. 
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TABLE I: Values of C a and Cf/K for various values of k in the case where the field 77 = (achieved by setting f3 to zero). 



There are three cases, depending on the value of k. The first case is n = k c . At this value, the potential energy 
is exactly zero for all separation. Indeed, the argument of both Bessel functions takes the same value. Furthermore, 
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the constants C a and Cf/n are equal, as can be see numerically from Table 1 and also analytically from the following 
argument. At n — n c , we can write the free energy in a form due to Bogomol'nyi psj : the field equations can then be 
written as first order differential equations: 

ip' - (™+KA)(p = and K^A' + ^j + 1(1 - <^ 2 ) = 0. (26) 

Linearizing the first equation, we obtain /' = —a; substituting Q into that equation, we find that the second 
Bogomol'nyi equation is satisfied if V%Cf — C a . This argument is only valid at n c because Il2(il) are not valid for other 
values of k. The potential is shown in Figure [3 

The second case is n > k c , corresponding to (conventional) type II SC. In this case, we note that Ko(y/2nd) falls off 
more rapidly than K${d) so the positive term in 125fl will dominate over the negative one for large enough separation, 
no matter what the values of the constants C a and Cf/n. Therefore, as is well known, in conventional type II SC 
vortices are repulsive and give rise to stable vortex lattice. The resulting potential is displayed in Figure [21 for n= 1. 

The third case is K < k c , corresponding to type I SC. Here, the situation is reversed: Ko(d) falls off more rapidly than 
Ko(y/2nd), and we conclude that the long-range potential is attractive. Thus, if we start with an initial configuration 
formed by a number of widely-separated vortices, ultimately they will collapse into a single vortex of large winding 
number. This case is illustrated in Figure[2]for k = 0.6 The result (|25[1 and the above discussion are in full agreement 
with |3, and pr| . 

Let us now turn our attention to the SO (5) case. The coefficients C a , Cf and C v are displayed for a variety of k 
and (3 in Table 2. 
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0.000 


3.0 


0.00 


0.362 


41.688 


0.000 




0.881 


1.419 


2.401 


0.049 




0.81 


0.362 


45.742 


0.083 




0.882 


1.422 


2.604 


0.081 




0.82 


0.364 


52.005 


0.119 




0.884 


1.428 


3.382 


0.122 




0.83 


0.364 


61.320 


0.145 




0.886 


1.434 


4.765 


0.152 




0.85 


0.368 


74.346 


0.185 




0.888 


1.440 


6.915 


0.177 




0.87 


0.371 


88.687 


0.216 




0.890 


1.446 


10.044 


0.198 




0.90 


0.379 


132.886 


0.249 




0.900 


1.481 


52.036 


0.281 













TABLE II: Values of C a , Cf/n, and C v /k for various values of k and j3 for m = 1 in the case of the SO(5) model. 



The possibility of an AF core gives rise to the surprising results mentioned at the beginning of this section. As has 
already been noted, if the core of the vortices is normal, there is no difference between the energetics of the SO (5) 
model and the SO(2) model, since the third term of i|2*5*|) is zero. Thus, the interaction energy of a pair SO(5) vortices 
(of winding number 1) is unchanged from that described above, to the left of curve A in Figure □ To the right of this 
curve, however, the r] field is nonzero and the third term makes an attractive contribution to the interaction force, 
either lessening the degree to which the vortex energetics behave in a type II manner, increasing the degree to which 
it behaves in a type I manner, or (the most interesting possibility) changing the vortices from repulsive to attractive. 

We can read at which point this occurs from (|25|l . The force between widely separated vortices depends on which 
of the fields is longest in range: the gauge field, which provides a repulsive force, or one of the scalar fields, whose 
forces are attractive. The range of the gauge field is 1 , while that of the rj field (the longest range of the two scalar 
fields, assuming it is nonzero in magnitude) is (VI — /3k) _1 . These forces are equal for a certain (3, which we denote 
(3*: 

/T = 1-1/k 2 . (27) 

The long-range vortex potential will be repulsive if (3 < (3* and attractive if (3 > j3* . 

Note that the caveat "assuming it is nonzero in magnitude" is essential: if the range of rj is longer than that of A, 
but Crf — 0, it will obviously not have an effect on the vortex interaction. 

A couple of case studies will help illustrate the situation; it is useful to see how the interaction energy varies with 
P for fixed k (see Figure|3J). In Figure^, the case k = 3.0 is displayed. (This corresponds to the second of the three 
cases mentioned near the beginning of this section.) The first two curves (/? = 0.0 and (3 = 0.8) coincide, since the 
core has yet to become AF (as can be ascertained from Figure ^ since these points lie to the left of Curve A). The 
subsequent curves show a decrease in the potential as the rj field turns on, along with its attractive contribution to 
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FIG. 3: Intervortex potential for k = 1 (a) and re = 3 (b) for m = 1, and for re = 1 (c) and re = 2 (d) for m = 2. 



the vortex energy. Eventually when /3 ~ 0.92, the potential goes from repulsive to attractive, in agreement with <|27|l . 
and also with 

In Figure Oa, the case k = 1.0 (the third case mentioned near the beginning of this section) is displayed. The first 
two curves, corresponding to parameters lying to the left of Curve A in Figure^ coincide, since the core of the vortices 
is normal; subsequent increases in [3 produce a drop in the interaction energy, as expected since r\ now contributes. 
It is important to note, however, that for re = 1.0 the superconductor is type I starting at (3 C = 0.33. We would thus 
expect the vortices to attract for (3 > 0.33. This is not the case: as can be seen from Figure[3Ji, the attraction begins 
only around f3 ~ 0.883, somewhat to the right of Curve A ((3af = 0.881). (That it is somewhat to the right of Curve 
A is simply due to the fact that r\ (at the center of the vortex, for example) evolves continuously with (3 and is strictly 
zero up until Curve A; it is therefore small and will have a negligible effect on the vortex interaction until it has time 
to grow to an appreciable value (see the values of C v given in Table 2).) Thus, there is a zone (between curves A and 
C) wherein the superconductor is type I, but wherein nonetheless vortices repel one another: this is the identity crisis 
referred to earlier. The point is that in this region the r\ field would produce an attractive long-range interaction, but 
it has not yet attained a nonzero value in the core and therefore has no effect on the vortex interaction. 

The situation is even stranger if we consider vortices of winding number 2, for which Curve B delineates normal 
(to the left) vs. AF (to the right) cores. For re = 2.0 (Figure^), the vortex turns AF at (3af — 0.562, but vortices 
remain repulsive until the value corresponding to the type I-type II transition {[3 C ~ 0.78). For re = 1.0 (Figure |3J;) , 
the vortices develop an AF core at roughly (3af — 0.659 but remain repulsive until just to the right of Curve C, as 
expected given that Curve C separates type II from type I superconductivity. However, at a lower value of re we would 
find that the interaction energy of 2-vortices in Region 6 is repulsive, which is incommensurate with the fact that the 
superconductor is type I there. 

Note that these "identity crises" only occur for certain values of the parameters, and in particular for re < 2.25, so 
this does not appear to be relevant to HTSC (supposing that the SO(5) model provides a good description of these 
systems), since all known HTSCs have quite large values of re (k > 50). However, in the case of superconducting 
cosmic strings and optical vortex solitons, there are no restrictions on the value of re and the effect may be relevant. 

It is interesting to note that the results obtained above complement those obtained in previous work where the 
conditions under which SO (5) superconductors are type I or type II were determined. In Ref. |22| . the distinction 
between type I and type II behaviour was determined by analysing the energy per unit flux of vortices as a function of 
their winding number. In Ref. |23| . this distinction was obtained by studying the energy density of a surface separating 
superconducting and normal regions at critical applied magnetic fields. The results stemming from both methods 
are in agreement and reveal, unexpectedly, that superconductors whose re values would normally be associated with 
type II superconductivity may in fact be type I superconductors. The surprising region is re > 1/V2 (conventionally 
type II) and to the right of Curve C in Figure ^ The underlying reason is the possibility of an AF vortex core, or 
equivalently, the fact that an applied magnetic field induces a transition between SC and AF (not normal) phases. 
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IV. CONCLUSION 



In this paper we have derived an expression for the potential energy of two widely separated vortices within the 
framework of a model with two order parameters, one of which is complex and attains a nonzero expectation value. 
Examples of physical contexts where such a situation occurs are bosonic superconducting cosmic strings, optical vortex 
solitons and the SO(5) model of HTSC. The approach used was a point vortex approximation, due to Speight [25| . 
Our starting point was a linearization of the free energy l|12ll to which source terms (I13f) were added; the form of the 
sources was given in H19I20|) . The main result is the two- vortex interaction energy, l|25|) . The constants appearing in 
this expression were chosen so that the long-range fields of a vortex in the point vortex approximation agree with the 
full nonlinear model. 

The behaviour of the vortices is for the most part in agreement with previous work, described in Section 2, namely, 
that for fixed k, type II behaviour (stable vortex lattice) turns into type I behaviour (unstable lattice) for (3 sufficiently 
close to 1, provided that the vortices have AF cores. 

This caveat, although apparently innocent, is actually quite important, and yields a surprising result. Namely, 
for certain regions in parameter space (regions 5 and 6 in Figure the superconductor is type I, yet vortices have 
normal (not AF) cores and repel one another, a behaviour normally associated with type II superconductors. It would 
appear that in such a case a vortex lattice would be metastable rather than stable (the free energy of a macroscopic 
non-SC region containing the same flux being lower than that of the vortex lattice). It is ironic and amusing that the 
surprising behaviour uncovered in [22, (type I SC for n > l/\/2) can be thought of as being due to the fact that 
vortices have AF cores, whereas the surprising behaviour uncovered here (vortices which repel one another in type I 
SC) is due to the fact that vortices have normal, not AF cores. 

While the attraction or repulsion of vortices is clearly of interest in any physical situation described by the above 
work, the unexpected behaviour uncovered here is unlikely to be seen in HTSC, as it occurs in the wrong region of 
parameter space. 

This work was supported by the Natural Science and Engineering Research Council of Canada and by the Fonds 
de Recherche sur la Nature et les technologies du Quebec. 
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